Semiclassical approach to discrete symmetries in quantum chaos 



We use semiclassical methods to evaluate the spectral two-point correlation function of quantum 
chaotic systems with discrete geometrical symmetries. The energy spectra of these systems can 
be divided into subspectra that are associated to irreducible representations of the corresponding 
symmetry group. We show that for (spinless) time reversal invariant systems the statistics inside 
these subspectra depend on the type of irreducible representation. For real representations the 
spectral statistics agree with those of the Gaussian Orthogonal Ensemble (GOE) of Random Matrix 
Theory (RMT), whereas complex representations correspond to the Gaussian Unitary Ensemble 
(GUE). For systems without time reversal invariance all subspectra show GUE statistics. There are 
no correlations between non-degenerate subspectra. Our techniques generalize recent developments 
in the semiclassical approach to quantum chaos allowing one to obtain full agreement with the two- 
point correlation function predicted by RMT, including oscillatory contributions. 



1 Introduction 

According to the random matrix conjecture [1] the energy spectra of quantum chaotic systems have 
universal statistical properties that depend only on the symmetries of the system and agree with 
predictions from random matrix theory (RMT). Usually one considers systems without geometrical 
symmetries. In this case one obtains agreement with the Gaussian Unitary Ensemble (GUE) if 
time reversal invariance is broken. The statistics of time reversal invariant systems agrees with the 
Gaussian Orthogonal Ensemble (GOE) if the time reversal operator squares to one, and with the 
Gaussian Symplectic Ensemble (GSE) if it squares to minus one; the latter situation can occur in 
spin systems. 

In this paper we want to consider individual chaotic systems with discrete geometrical symme- 
tries. The energy levels of these systems fall into subspectra associated to irreducible representations 
of the underlying symmetry group. For example, in systems with a reflection symmetry the sub- 
spectra correspond to eigenfunctions that are either even or odd under reflection. Moreover, if the 
system is chaotic each subspectrum individually obeys RMT statistics. One remarkable property, 
first observed by Leyvraz, Schmit and Seligman [2], of symmetric systems is the ability to generate 
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GUE statistics within certain subspectra regardless of the time-reversal properties of the fuU sys- 
tem. A partial explanation of these observations was given by Keating and Robbins [31 H] who used 
semiclassical analysis to evaluate (via its Fourier transform) the two-point correlation functions 

associated to each subspectra. Here p{E) = ^„ 6{E — En) is the corresponding level density, p is 
the mean level density, and (•) denotes an energy average. The GUE and GOE predictions for R{e) 
are [5l [6] 
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Here we have split the correlation function into oscillatory and non-oscillatory contributions. The 
GUE result is proportional to \ whereas the GOE result leads to an infinite asymptotic power 
series in ^. The terms involving odd powers of i are meaningful if e is taken with a small positive 
imaginary part. Keating and Robbins used the so-called diagonal approximation to recover the 
leading non-oscillatory term proportional to ^ (corresponding to the linear term in the Fourier 
transform). Depending on the time reversal properties of the system as well as the type of the 
representation they obtained as predicted by the GUE or — ^ as predicted by the GOE. 

In the present paper we show this agreement with RMT carries over to the full correlation 
function by building on recent results for systems without geometrical symmetries. The starting 
point for all such work is the Gutzwiller trace formula [7] 

p{E) ~ -p{E) + ^Re V TaFae'^'^^'^y^ (4) 

a 

which relates the quantum density of states to a sum over the actions Sa{E) of the periodic orbits 
of the corresponding classical system, weighted by the primitive period Ta and complex stability 
amplitude Fa (which contains the Maslov index). Convergence can be ensured by taking the energy 
E with a large enough imaginary part which limits the contributions from long periodic orbits. By 
directly inserting the trace formula into ([T]) one sees that spectral correlations are determined by 
correlations between periodic orbits [8] . The evaluation of these correlations allows one to obtain 
all non-oscillatory contributions to R{e). The first term is given by the aforementioned diagonal 
approximation [9l [10] taking into account pairs of identical or (in time-reversal invariant systems) 
mutually time-reversed orbits. The remaining terms arise from pairs of periodic orbits that differ by 
their connections inside close 'encounters' |ll ^ ll2 ^ [T3 ^ ll4 ^ ll5j. For example, one orbit could contain 
a self crossing with a small intersection angle while the other orbit narrowly avoids the crossing. By 
analysing the nature of such encounters one can also show why for broken time-reversal invariant 
systems there are no further non-oscillatory terms beyond the diagonal approximation. 

Currently our knowledge of orbit correlations is not sufficient to obtain the oscillatory com- 
ponents directly, instead one must take an indirect approach by relating long and short orbits 
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through a 'resummation' or 'bootstrapping' procedure. This technique can be apphed directly to 
the periodic orbit sum Q [16j , however it turned out more convenient to consider the resummation 
of pseudo-orbits (collections of periodic orbits) instead. The resulting Riemann-Siegel lookalike 
formula |17l \TE[ [T9] manifestly preserves the unitarity of the quantum time evolution and was com- 
bined with a so called generating function [20] in order to access the oscillatory contributions to 
the two-point correlation function [21j . Inspired by field-theoretic approaches the authors of P2] 
were then able to relate contributions from sets of pseudo orbits with small action differences to 
terms arising in the nonlinear sigma model. This resulted in the reproduction of the full GOE and 
GUE expansions for systems with and without time reversal invariance. 

We will combine this approach with results from representation theory to investigate systems 
with discrete geometrical symmetries. Our techniques have interesting analogies to the work of 
Boris Gutkin which recovers the non-oscillatory contributions to R{e) for "cellular" billiards [23]. 
Correlations between different subspectra are also considered in [23j. 

The paper is organized as follows. In Section[2]we will outline some of the basic principles behind 
the application of representation theory in quantum mechanics. In Section[3]we will then introduce a 
generating function determining the correlations inside each subspectrum. This generating function 
and the corresponding correlation function will be evaluated in Sections |4] and [5] using correlations 
between orbits as well as results from representation theory. In Section |6] we will generalize to 
cross-correlation functions between different subspectra as well as the full correlation function of a 
symmetric quantum system. Conclusions will be given in Section [7j 

2 Symmetry in quantum mechanics 

To get started, we give a brief outline of parts of representation theory that are helpful for under- 
standing symmetries from a semiclassical point of view. For details see e.g. [25l l26| [271 [28] . 

Classical symmetries are characterized by a group G of symmetry operations that leave the 
system invariant. For example if a classical Hamiltonian H{q,p) is invariant under rotations by 
27r/3 the symmetry group is given by G = C3 = {e,g,g^}, where g denotes rotation by ^ and e is 
the identity. Figure [l] shows a billiard system with this symmetry. 




Figure 1: A billiard system that is invariant under rotations by 27r/3. 

If the group G consists of point transformations then in quantum mechanics it induces trans- 
formations of wave functions through the definition U{g)ip{r) = ip{g~^{r)), g & G. The unitary 
operators U{g) form a representation of the group G because they satisfy U{g2) U{gi) = U{g2gi) 
for all gi,52 G G. In a symmetric system the transformations U{g) commute with the Hamiltonian, 
\U{g)^H] = for all g £ G, and G is said to be a symmetry group of the Hamiltonian. One 
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can consider also more general symmetry groups for which the operators U{g) are not related to 
coordinate transformations (for example, the Hecht Hamiltonian in [3] has a symmetry in angular 
momentum space). 

In a quantum system with a discrete symmetry group one can split the spectrum into appropriate 
subspectra. For example, if the symmetry group is C3 then one can define three subspaces by 
requiring that the wave functions are periodic with respect to rotations by 27r/3, up to a phase 
factor whose third power is one: ^{r, (p — 2tt/3) = e^'^^"'^^'ip{r, (p) where q = 0, 1, 2. More formally, 
it follows from [U{g),H] = that one can find a basis in which the action of U{g) on a particular 
basis function is determined by the irreducible unitary matrix representation M^^^g) (see |25j-|28j) 
as follows 

Uig)\a,i,n) = ^Mjf{g)\a,j,n), H \a,i,n) = E^^^a,i,n) . (5) 
i=i 

Here a labels the different irreducible representations, i the components within this representation, 
n labels basis functions with the same a and i, and Sa denotes the dimension of the representation. 
Hence when U {g) acts on the entire Hilbert space it is block-diagonal, with each block an Sa x Sa 
matrix given by M^'^\g). Moreover the Sa eigenfunctions that correspond to this block all have 
the same energy and so the eigenvalue E^"^ is Sa-fold degenerate. One can hence decompose 
the spectrum of a symmetric system into subspectra that correspond to the different irreducible 
representations of the symmetry group. 

Note that the irreducible representations of a group satisfy Yla^'a — \G\ where |G| is the 
number of elements in G. In the example above the symmetry group C3 has three one-dimensional 
representations, labelled by a = 0,1,2, and given by the one-dimensional matrices M^'^\g^) = 
^2mNa/3, £qj, = 0, 1, 2. Hence the eigenfunctions of the Hamiltonian satisfy 

U{g'')i;'C\r, cP) = i;t\r, <P - 27riV/3) = e^-^^/^V^l") (r, (6) 

in agreement with the more informal approach above. For completeness we note that the projection 
operator onto the subspace a can be explicitly given in the form [251 [26] 

^" = ^E^^(5)trM(")(5), (7) 

satisfying Pq, |a', z, n) = SacH \a.\i,n) and = Pa- 

In the following it will be important that one can distinguish between three different types 
of representations. An irreducible representation is real if the matrices M^"\g) are real, or can 
be made real by a simultaneous similarity transformatiorj^ The representation is pseudo-real if 
it is not real but all matrices M^"\g) have real traces. In the remaining case where not all 
matrices have real traces the representation is complex. Complex representations come in complex 
conjugate pairs, because for every complex representation M^°^\g) there exists a complex conjugate 
representation M^^\g) = M^"\g)* , which is inequivalent to M^^^g). If a quantum system is 
invariant under complex conjugation, which corresponds to the time reversal operation for non- 
spin systems, this leads to an additional degeneracy between the spectra of complex conjugate 
representations. For example, for a time-reversal invariant system with symmetry group C3 one can 

^Irreducible representations that are related by a simultaneous similarity transformation M^"\g) — SM^"'^ {g)S~^ , 
g £ G, are equivalent. 
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use complex conjugation to obtain from every solution that satisfies \lj{r, <f> — 27r/3) = e^'^^^^ip{r, (f)) 
another solution that satisfies il^{r, (p — 27r/3)* = e~^'^^^^'tp{r, (/>)* = e^'^*/^'i/'(^) 4')* has the same 
energy. For real and pseudo-real representations the complex conjugate representation is equivalent 
to the original on^ 

An important property of the representation matrices M^'^\g) that we will use in the following 
is the group orthogonality relation that follows from Schur's lemmas. It states that 



^E^5(^X(^r = ^^^^- (8) 



In the semiclassical approach the representation matrices usually contribute through their traces 
Xa{g) = tr M(")((7), called characters. By taking traces in ^ one obtains the orthogonality relation 
for characters ^ 

VQ\^Xa{g)X*p{9) = ^aH- (9) 
' ' 9&G 

We will also need the following identities for group averages of characters 

7777 y'xa (05^5) = —Xa{ab'^), (10) 

^,^Xa{ag)Xa{bg''^) = —Xa{ab), (11) 

I I g&G 

77^ y]xQ(a5^5""^) = —Xa{a)Xa{b). (12) 
1*^1 gee '^'^ 

Here o and b are arbitrary group elements. In the first identity we respectively have Cq = 1,0,-1 
for real, complex, and pseudo-real representations. The second and third identity hold regardless 
of the type of representation. These identities are proved in the appendix of [29]. They can be 
obtained by manipulating a simpler identity from [26j and the group orthogonality relation ([s]). 

In this article we will show that for spinless time reversal invariant systems the subspaces 
associated to real representations obey GOE statistics, and subspaces associated to complex rep- 
resentations obey GUE statistics in spite of time reversal invariance. Pseudo-real representations 
lead to GSE statistics |30j but will be excluded in the present paper. For systems without time 
reversal invariance all subspaces obey GUE statistics. Different subspectra are uncorrelated apart 
from the degeneracies mentioned above that are related to time-reversal invariance. 



3 Generating Function 

We are interested in the statistics of the energy levels in a subspace a, where the Sa-fold 
degeneracy has been removed. The subsequent level density Pa{E) = X]„<5(-E — -E^"^) can be 
accessed from the following trace 

Pa(E) = lim Imy^ = lim Im f — trP„— —] , (13) 

' TT E+^E ^E+-Et^ TT E+^E \Sa E+ - H J ' ^ ' 



^ An additional two-fold degeneracy is also present in pseudo-real representations if the quantum system is invariant 
under complex conjugation due to a different mechanism [271 128j . 
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where Pa (see Eq. ([7|) is the projector onto the subspace a and = E + ir] has a small positive 
imaginary part that is taken to zero. Similarly the correlation function 



Raie) 



1 



Pa[E + 



Pa{E 



1 



(14) 



can be obtained as the limit for e"^ 

1 / 1 



e + zry — )• e of the real part of the complex correlator 



Ca(e+) 



-trP„ 

•5q E + 



2-Kpa 



H Sc 



' E- 



H 



(15) 



Here pa is the mean level density of the subspectrum. It is related to that of the full spectrum by 
Pa ~ Sap/\G\ [31r The random matrix predictions for the complex correlation function are given 
by Eqs. ^ and ([3), in each case dropping the restriction to the real part. 

To access the level density using Riemann-Siegel resummation we must first relate the trace of 
the symmetry projected resolvent to the corresponding spectral determinant Aq,(£') that vanishes 

at the energies eI°'\ The connectfon is given by the usual relation that the trace of the resolvent 
is the logarithmic derivative of the spectral determinant. Hence we hav^ 




and vice versa 



/S.a{E^) oc exp 



1 1 



-tr 



Pa 



E' -H 



dE' 



d AaiE) 



dE' Aa{E') 



E'=E 



If we use this idea for both traces in (15) we can then access the complex correlator by 

2- 



Ca(e+) 



I) 



deAdes 

where the generating function Za is defined by 

' Aa{E + ec/2Tipa)Aa{E - 60/27: Pa) 



(16) 



(17) 



(18) 



Za{^A, e_B, ec, ^d) 



Aa{E + eA/27rpa)AaiE - eB/2TTpa) 



(19) 



and (II) denotes the matching conditions €a, ^B, = e^- 

We now want to derive a semiclassical approximation for the spectral determinant and thus 
for Za- For this purpose we need to find a semiclassical expression for the trace of the symmetry 
projected Green's function. Fortunately this has already been obtained by Robbins [3j (see also 
Seligmann and Weidenmiiller |i33]) and is given by 



1 



-tr 



Pa 



1 



E+ -H 



iSa{E+)/h 



(20) 



*Here we have translated the result of [31] to our notation. In contrast to [31] our pa contains each degenerate 
level only once. 

^In general, the trace of the resolvent and the determinant need to be regularised, see [35] ■ This leads to a 
regularisation dependent prefactor of the determinant that it is not relevant for the following semiclassical calculations. 
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Here ga{E~^) is the smooth part of the trace with the l/s^ prefactor. It is related to the mean 
density of states by Im^o- (£'"'") = —Trpc,{E~^). Expression (20) contains the same primitive orbit 
period Ta, action Sa and stabihty amphtude Fa as the usual Gutzwiller formula Q, however now 
the sum is over all orbits a which are periodic in the fundamental domair^ but not necessarily 
periodic when unfolded to the full system. However their final point is always related to the initial 
point via one of the symmetry operations ga- The weight then contains the character Xa{ga) that 
corresponds to this symmetry operation in the representation a. 

This may be illustrated using the stadium billiard which has a four-fold symmetry given by 
the group 82(^82 = {e,rx,ry,rxry}, generated by the two refiection operators and Vy. The 
fundamental domain consists of a quarter of the stadium with specular reflection conditions at its 
boundary (shown in grey in Figure [2]) . Periodic orbits in the fundamental domain may be unfolded 
to the full system as shown in Figure [2j In this example, a periodic orbit that retraces itself in 
the fundamental domain where it strikes both the x and the y axis (Figure [2]^a)) will traverse 
both these axes in the full system and finish in the domain associated with the r^Vy group element 
(Figure [2|^b)). Hence the corresponding weight is Xa{fxTy). We may view the group generators as 
symbols, allowing each section of an orbit to be assigned a particular symbol sequence and therefore 
a particular group element comprised of these generators. 





Figure 2: (a) A self-retracing periodic orbit in the fundamental domain is unfolded to (b) a partial 
orbit in the full stadium. 



If we insert (20) into (16) we obtain the following semiclassical expression for our spectral 
determinanll3 



A, 



\ a. 



,iSa(E+)/h\j 



(21) 



where Na{E) is the mean counting function obtained by integrating pa{E). The sum is now over 
pseudo-orbits (finite collections of periodic orbits) A with cumulative actions Sa and amplitude 

^ We omit the contribution from orbits confined to tlie boundary of the fundamental domain. These orbits involve 
multiple repetitions of a primitive periodic orbit. Due to the exponential proliferation of orbits with increasing period 
their contribution is negligible. 

'^We neglect in the following the contributions of non-primitive periodic orbits to the trace formula since they are 
exponentially suppressed in the limit of long periods. Furthermore, the real part of Qa is regularisation dependent 
and included in the proportionality factor. 
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factors G^^\ given as follows 



G 



(a) 



n 



Sa — ^ ^ T^aSa 1 



riA 



(22) 



Each Ua runs from to oo, but ua is always finite (and can be zero). For energies with negative 
imaginary parts one uses the relation Aa{E~) = [Aa{E~^)]* . All information about our irreducible 

subspace is contained within the factor G^T^ and more specifically the characters within this factor. 



We also note the similarity between the expression (21 ) and the symmetry reduced dynamical zeta 
functions defined in 



At this point our motivation for employing the spectral determinant becomes clear. As it 



stands, the expression (21) suffers from exactly the same convergence issues as the Gutzwiller 



formula. However Aa{E) obeys a functional equation which, via analytic continuation, can be used 
to show that truncating the sum in (|2l|) at half the (rescaled) Heisenberg time ^ — ^irhpa gives 
a contribution to Aa{E) approximately equal to that of the complex conjugate of the remaining 
orbits [19]. In our context this resummation procedure, known as the 'Riemann-Siegel lookalike' 
formula [n\ [T8] leads to an improved semiclassical approximation of the form 



A„(£) 



E ( 



-1 



G 



exp 1-Sa{E) - inN^iE) \ + 



c.c. 



(23) 



where Ta = dSA/dE is the cumulative period of pseudo-orbits in ^. In the semiclassical limit 
— )• oo and one recovers all contributions. We now follow [_22j by inserting (23) into the spectral 



determinants in the numerator of our generating function (19). This gives rise to four terms 
with oscillatory components in the numerators of the form exp{ziziTr[Na{E + ec/^npa) — Na{E — 
ez)/27rpa)]) and e-x.p{±.m[N ^{E + ec/27r/)a) -|- Na{E — e£)/27r /)„)]). However those with additive 
phases oscillate rapidly and are assumed to give zero contribution after averaging over the energy 
E. Thus our semiclassical generating function may be written as the sum of two terms 



Za{eA,eB,ec,eD) = Z^\eA, es, ec, en) + Z^\eA,eB,e*c,e 



D) 



(24) 



with the first given by 



\ A 



B 



X e 



-iiTNc{E+ec/2-Kpc,) 



c 



X e 



i-KNa{E-to/2lTpa) 



TlD M'^)* „-iSD{E-tu/2-Kpa)/h 



Y.{-lT-G^S^*e 



(25) 



D 



and the second obtained by taking the complex conjugate of the third and fourth line, or equivalently 

Z(2)(e^, es, e*c, = zj^\eA, es, -eh -e*c) . (26) 
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To simplify our expression we expand the mean counting function and action to obtain 



Na{E ± e/2^/>„) ~ No,{E) ± e/27r (27) 

and 

S{E ± e/2Tip^)/h ~ S{E)/h ± T{E)e/T^^^ (28) 
which are then inserted into the generating function to give the fohowing expression 

^(1) ^ gj{eA+es-ec-fD)/2 / ^ (-1)""^+"^ 

^ gi(SA-S's+5c-Si3)/rtgj(TAeA+Tses+Tcec+T'DeD)/T^"'y ("29) 

Contributions to the generating function will occur when the action difference /S.S = Sa + Sc — Sb — 
Sd is at most of the order of h; summands for which the action difference is large in comparison to 
h are expected to wash out after averaging over the energy. Systematic contributions will obviously 
arise when the same orbits (modulo time-reversal) occur in both AuC and BUD, extending Berry's 
original 'diagonal approximation'. The remaining 'off-diagonal' contributions will arise if the orbits 
in B L) D differ from those in j4 U C only due to their connections inside encounters. Using these 
ideas the generating function of systems without symmetries was evaluated in \20\ [211 122j . The 
following sections extend this approach to systems with discrete symmetries by evaluating what 
effects the characters Xaida), associated with each periodic orbit, have on these contributions. 

4 Diagonal Approximation 

We start by evaluating the contribution from those orbits which have the same action in time- 



reversal invariant systems. For this it helps to rewrite (29) by converting each of the four pseudo- 



orbit sums back into an exponentiated sum over periodic orbits, i.e. 



ZW = e^(e.+e.-ec-ez,)/2 ^,^p FaXa{9a)e''^/' flc^ exp F,* xMe-"''^' fno^ ^ (30) 



with 

The exponentials can then be re-expanded in terms of Taylor series to give 



j CO ^ n m 



bi,...,hm. 



(32) 

To implement the diagonal approximation we assume that all contributions may be neglected unless 
the actions satisfy Yli = Xlj '^bj ■ For this to occur a periodic orbit Oj must be matched to a 
periodic orbit bj that is either identical to Oj or mutually time reversecFl hence we can set n = m. 



Again we neglect orbits that are just repetitions of shorter orbits. Furthermore, it is sufficient to consider only 
the case where all ai are different, because the other cases lead to negligible contributions. 
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Given that there are n! ways of matching the Oj to bj we can then cancel one of the factorials in 
the denominator. Furthermore, if the orbits are identical we can replace the character Xaidbj) by 
Xaidai): whereas if they are mutually time reversed the corresponding group elements are mutually 
inverse and the characters are mutually adjoint (due to the unitarity of the representation). This 
allows us to replace Xa iObj ) by Xa {da^ ) — Xa {Oa^ ) ■ The resulting expression can then be written in 
the form 



Z 



(1) 

a,diag 



P«(eA+es-ec-eD)/2 



LXJ _j It 



(33) 



with 



l^a{a) = Xa{9a)[X*a{9a) + XaiQa)]; (34) 

Note that KQ,(a) replaces the factor 2 for time-reversal invariant systems without symmetries. Rear- 



ranging the sum and product in ( 33 ) we recover the standard Taylor expansion in n of an exponential 
yielding 

.(1) 



Z 



a,diag 



\ a 



\r? \2 fa fa 
\^a\ JaCJBD 



(35) 



Now following the arguments in [1] we assume that the stability amplitudes Fa and the group 
elements ga of the orbits are uncorrelated. This allows Ko(a) to be separated from the orbit sum 
and averaged independently over the group. Due to ergodicity this group average is uniform as the 
probability of finding the end points of an unfolded orbit in any particular copy of the fundamental 
domain becomes uniform as T — )• oo. Therefore we may write 



.(1) 



■ jdiae 



(36) 



where 



^diag 



{Xai9)[X*ai9) + Xaig)]) , 



1 



E 

9 



Xai9)\'^ + 



1 

W\ 



9 



Xaigf 



(37) 



and the nature of the exponential allows us to change the factor of Kq'"^ to a power. In (37) the first 



summand is simply 1 due to the orthogonality relation ^ whereas the second summand depends 
on the type of representation. If a is real then the characters are real and we obtain the same result 
as for the first summand. However if a is complex then the complex conjugate of a character is the 
character of the complex conjugate representation and the second summand vanishes due to 
Hence 

' I " ^'""^ , (38) 
1 a complex. 



,.diag 



Now, by using the Hannay and Ozorio de Almeida [9j sum rule we can approximate the sum over 



periodic orbits in (36), weighted by stability amplitudes \Fa\^ , by an integral of the form J"^ 



Here Tq is the minimum period from which the orbits tend to behave ergodically and after scaling 
with the Heisenberg time takes the lower limit of the integral Tq/T^^ — )■ in the semiclassical limit. 
Computing this integral (see [22j for details) then leads to a semiclassical generating function 



Z^^^ = e 

a,diag 



«(eyl+<:s-<:c-eD)/2 



(ec + eB){eA + ep] 
i^A + eB)(ec + es) 



(39) 
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(2) J l 

where ^^^^ can be obtained from the relation (26). This aUows us to obtain our complex correlator 
through differentiating our generating function as follows 



r (r\ - o l^W ,7(2) 

•-^a.diaglt; — Q^^Q^^ [^a,diag "i" ^a.diag 



(ii: 



1 

2' 



Finally, inserting ( 39 ) into the above relation when a is complex we achieve 



a,diag 



'2e2 + 2e2 



(40) 



(41) 



which is exactly the diagonal approximation predicted by RMT for the complex correlator in the 
GUE case Q. It is worth noting that this happens precisely because the second term in (37), 



attributed to the correlation between a periodic orbit and its time-reversed partner, averages to 
zero. Hence, as remarked in [3] we observe in complex representations a mechanism akin to the 
introduction of an Aharonov-Bohm flux, in which time-reversed periodic orbits exist but their 
attributed phase factors average to zero [36]. In contrast, for real subspaces the coefficient k'^'^^ = 2 
means 

Ca,diag(e) = -\ (42) 



which corresponds to the diagonal term in the GOE complex correlation function. 

Conversely, in systems where the classical time-reversal symmetry is broken those time-reversed 



periodic orbits are no longer available and so the coefficient (37) becomes 



= {Xa{9)x*a{9)), 



a real 

a complex. 



(43) 



and we obtain the GUE complex correlation function in both real and complex subspaces. In 
summary, with the diagonal approximation we showed agreement with the RMT predictions up to 
leading order e"^^ and it remains for us to validate all other orders. 



5 Off-Diagonal Contributions 

Further contributions to Za^ with small action differences arise if the orbits m BUD closely follow 
those in ^ U C except in so-called encounters where several orbit stretches come close (see Figure 
[3]). By changing the connections inside these encounters (which may lead to orbits splitting or 
merging), we can relate the orbits m. AuC to those in BUD. The possibility to switch connections 
in this way is a consequence of hyperbolicity. In contrast the 'links', i.e. the long parts of the orbits 
connecting the encounter stretches, are almost the same in C and BU D modulo time reversal. 



We have to evaluate the contribution Z^^^^ to the generating function ( 29 ) arising from AU C 



and B U D differing by their connections in encounters, i.e., 

^ gj(€A+es-ec-ei5)/2 ^ G^?"* G^^* (- 1)"'^+"° 

A,B,C,D 
diff. in enc. 

^gjA5/Sgi(r4eA+Tses+Tcec+TD€D)/T^"'\ ^ (44) 
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Figure 3: (a) A topological configuration space view of pair of periodic orbits with a single en- 
counter, (b) Collections of bunched periodic orbits forming pseudo-orbits differing in multiple 
encounters (from [37]). 



where A^* = Sa + Sc — Sb — Sd- To account for the possibility that some orbits in B L) D are 
obtained by changing connections in AU C whereas others are just identical Z^l^^ still has to be 

multiplied with Z^\ii^^- Altogether we thus obtain the following expression for Za^ 

4^^=4iaJl + ^i:iff) • (45) 



a,diag \ ' a, off J 

At present we only consider systems with time-reversal invariance. Systems without time- reversal 



invariance will be discussed in section 



5.3 



The pseudo-orbit quadruplets that contribute to Z^^^g 



can have many different topologies. Following [22] (to which we refer for details) these topologies 
may be classified by using the notion of 'structures' of quadruplets {A, B, C, D). These structures 
are characterized by the number of encounters, the number of orbit stretches participating in each 
encounter, and the ordering of these stretches along the orbits. In systems with geometric sym- 
metries, the contribution from each structure may be evaluated similarly as for the non-symmetric 
systems in [22] provided we account for the group characters associated to each periodic orbit 
and relevant factors of |G| and Sq that arise from considering orbits in the fundamental domain. 
We discuss in the following the modifications of the calculations in [22] due to the presence of 
symmetries. 

First we separate each pseudo-orbit oi AuC into the constituent periodic orbits pi, . . . ,pn and 
each pseudo-orbit oi B U D into constituent periodic orbits qi, . . . ,qm- The combined amplitude 
factors from orbits in AU C can be considered approximately equal to those in B L) D. We then 
determine the action difference between both sets of periodic orbits which depends on the separation 
between the stretches involved in the encounters, with components s and u pointing in the stable 
and unstable directions in phase space. We thus require a probability density ■wti,T2,...t„('S, u) for 
finding in given periodic orbits {pi to Pn) V encounters with a given structure and given s and u. For 
systems without geometrical symmetries this density was evaluated using ergodicity and was shown 
to be proportional to where Q is the volume of the energy shell and L = Yla=iK'^) is 

the total number of encounter stretches (equivalently links), obtained from the number of stretches 
l{cr) in each encounter a. For systems with geometrical symmetries we must replace with the 
volume of the fundamental domain which leads to a factor of in comparison to the 
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density in [22]. Summation over structures and integration over s and u therefore yields 

struct ' I l(T V / p-^^..,^p^ 

/'''"■'''^""- .(.").-'">'V-"-..™»./<'>. ,4.) 

Here the division by 2^F!]^^i((T) avoids overcounting due to choices of {A, B,C, D) that can 
be described in terms of several equivalent structures and the coefficient Ka{pi, ■ ■ ■ ,qm) simply 
assembles the characters associated to each periodic orbit, i.e. 

Ka{Pl,---,qm) = Xa{gpi)--- XaigpJxUdqi) ■ ■ ■ XliOgJ- (47) 

As in the previous section the crucial step then is to replace this coefficient by an average over all 
possible group elements for each periodic orbit that are consistent with a considered structure. 

C"" = iXaigp,) ■ . . Xa{9pJXai9<li) ■ ■ ■ X*a{9<lJ) • (48) 

We will discuss later in more detail how this average is performed for a given structure. Note that 
if one replaces a periodic orbit by its time-reversed equivalent then this corresponds to a different 
structure in the present formulation. 

In summary, we have made the following changes compared to the non-symmetric case. In addi- 
tion to the factor k^'"''' and the factor of arising from the probability density wti,T2,...t„{s, u) 



we have also replaced Th by in the exponent of ( 46 ) . This amounts to a multiplication of the 
exponent by Th/T^^ = p/pa = \G\/sa- This can be taken into account by replacing all the ej by 
€j\G\/sa in the final result for the off-diagonal contribution to the generating function. 
In the non-symmetric case the final result was [22j 

7(1) _ V 7(1) ^T^p.p 7(1) _ (-1)"^+"^ rionc ij^AorC + egorP) . . 

^off ~ ^off.stmct' wneie '^off, struct ~ oVt/I TT -La \' ^ ^ 



struct 



Here each link that is included in an original orbit pi and a partner orbit qj gives rise to a factor 
with subscripts depending on whether pi belongs to A or C and whether qj belongs to -B or D. 
The subscripts in the encounter factor depend on the orbits containing the beginning of a specified 
encounter stretch. 

According to the discussion above we have the following modification in the symmetric case 

7(1) _ structi^iL-y ^7(1) _ struct„ L-V yW (r,n\ 

^a,oS- I'-'l \~ J ^off, struct - Z^^a ^off, struct VO'JJ 

struct ^ ^ struct 



where ^off struct given by (|49|) and is symmetry independent. 



As ■^off struct known we only have to deal with k^'^'^'^^. We will see that for time-reversal 
invariant systems k^^"^* makes sure that structures that actually require time reversal invariance 
are omitted if the representation is complex. In all other situations k^^^^'^^ cancels the factor Sa^~^ ■ 

if Q real, or a complex and structure does not require time-reversal invariance, 
if a complex and structure does require time-reversal invariance. 



struct I 



(51) 
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5.1 Real Irreducible Representations 



We start by evaluating kJ^''"''* for real irreducible representations in time-reversal invariant systems. 
To illustrate our approach we consider the orbit pairs introduced in and depicted in Figure |4j 





Figure 4: An orbit with two antiparallel encounter stretches, and its partner orbit obtained by 
changing connections inside the encounter. 

Here AuC only contains one periodic orbit p where two encounter stretches are almost antipar- 
allel. In the partner orbit q, included in B U D, the connections inside the encounter are changed 
such that one of the links, say the first, is reversed in time. Importantly, the group elements asso- 
ciated to p and q can be decomposed into the group elements associated to each link and encounter 
stretch. These elements indicate the symmetry operation that relates the copy of the fundamental 
domain in which the unfolded stretch or link starts to the copy in which it ends. If a stretch or link 
is changed slightly the group element stays the same and if it is reverted in time the group element 
is inverted. Thus in our example p has a group element gp = l2elie~^ (with e and li denoting 
encounter and link elements respectively whereas in q the first link is reverted in time such that 
9q = h^li^e^^. The corresponding k^''"''' is therefore 

C"" = ^ E xUl2elie-')x*aihdi'e-'). (52) 

However the encounter elements can be removed as we can decompose e as e = urj and use the 
invariance of the group average under multiplication by an arbitrary group element to substitute 
rjliT]^^ —7- li and i^^^l2i^ — )• h- This leads to 



Kr" = 42 E xa{i2h)xuyi') (53) 



\G\' 



which in the following sections we will show to yield 

In general for a particular structure k^''"''* may be written as 

= iXaigp,) ■ ■■Xa{9pJX*a{9q^) ■ ■ ■ X*a{9qJ) e,,...,ev ,h,-,lL ' (^^^ 

Here the group elements of the orbits gp-^, . . . , gp^ , (7^^ , . . . , gq^ are alternating sequences of group 
elements ei, . . . , ey associated to the encounters (or their inverses) and group elements . . . ,Il 
associated to the links (or their inverses). However, as in the example above we may split the 
encounters into the form Cj = I'i'qi and redefine the link elements using the invariance of group 



sums under multiplication. Thus all encounter elements may be dropped from (54). 

We now proceed to evaluate kJJ''"''*. Bolte and Harrison have encountered a related situation 
\38\ [29] on quantum graphs, where the group elements refer to the precession of spin around a 



periodic orbit, rather than unfolding periodic orbits in a system with a discrete symmetry. See also 



must not be confused with the identity element! 
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\39\ [13] for the case of flows. In these instances group averages were performed for pairs of periodic 
orbits which determine the non-oscihatory contributions to the correlation function, without using 
the pseudo-orbit approach for the generating function. Here we extend these methods to discrete 
symmetries and structures involving arbitrarily many orbits. 




BuD AuC 



Figure 5: For an encounter with three stretches the connections can be changed in two steps each 
affecting only two stretches. 

We will use the fact that to fully change the connections in an encounter of / stretches we 
need to perform / — 1 steps that just interchange the connections between two encounter stretches. 
Therefore to fully change connections in all V encounters and change the orbits in i? U D to those 
contained in ^ U C we need a total of Yl^=iiK'^) — I) = L — V such steps. For example, Figure 
4 depicts a structure involving 3 — 1 = 2 steps. There are different types of steps, depending on 
whether the stretches involved are parallel or antiparallel and whether they belong to the same 
orbit or different orbits. We will show that regardless of these options each step gives a factor of 
Sa~^, multiplication then leads to Sa^^ 



5.1.1 Anti-parallel encounter permutations 




Figure 6: (a) Black/full line: an orbit p £ AU C for which an encounter stretch (big arrow) that 
connects two links li and I2 has been highlighted. P denotes the remaining part of the orbit. 
Red/dashed line: an orbit q £ B L) D that traverses the link li in the opposite direction and I2 in 
the same direction as p. (b) Black/full line: p again. Red/dashed line: the orbit q' where one of 
the reconnection steps leading from q to p has been performed. 



To derive the factor we compare the quantity k^''"''' in Eq. (54) to the one for a different 
structure where one of the L — V steps leading from B L) D to Au C has been performed. Let us 
begin by selecting an individual encounter stretch of an orbit p £ AU C . This stretch connects two 
links whose group elements will be denoted by and I2 (and its group element is absorbed in li 
and I2 as discussed above) . Hence the group element associated to p will be of the form 

gp = Phh , (55) 

where the product of group elements has to be read from right to left. Here P is associated to the 
remaining parts of p. Its precise form is unimportant for our calculation. The links li and I2 also 
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show up in i? U D. We first consider the case where one of these hnks, say li, has been reverted in 
time and the two hnks form part of the same orbit q G B U D. This situation is depicted in Fig. |6j 
The other cases are treated later. If we denote the group elements of the remaining parts of q by 
Qi and Q2 then q has a group element of the form 

= {Q2h)[Qi\^) ■ (56) 

As illustrated by Fig. [6j the orbit p contains an encounter stretch following li and preceding I2 
(indicated by the big arrow). The orbit q has an almost parallel encounter stretch preceding the 
link and an almost antiparallel encounter stretch preceding the reverted link (the two crossing 
big arrows). In the decomposition of gq in (56) the locations of these stretches are indicated by 
brackets. 

Now we transform the structure to make BU D more similar to AuC, by changing connections 
between the two antiparallel encounter stretches singled out above. This replaces q by an orbit q' 
where one of the brackets in gg, say {Qi^l^^), has been inverted, leading to 

9,' = {Q2h){hQi) ■ (57) 

In q' the links with group elements li and I2 follow each other directly, as in p. The structure in 
which q' replaces q is hence simpler than the original structure and now one less reconnection step 
is needed to go from BiJ D io Au C . 

We now want to relate the corresponding coefficients k^''"'''. By inserting the orbit decomposi- 
tions (55) and (56) for the original structure into ( |54[ ) its coefficient can be represented as 

C"" = (Xa(fe)Xa(fe)...Xa(^'^2/i)x;(5,Jx;(5,J...X;((Q2/2)(Qr'^r'))> (58) 
with averages running over all link elements. In particular /t^''"''' involves the sum 

E Xa{Pl2h)xl{{Q2l2){Q^H^^)). (59) 

We now want to bring this sum to a form that contains gqi instead of gq. For this purpose we use 
again that group averages are invariant under multiplication by an arbitrary member of G. This 
allows us to replace h — >• x~^li and I2 —5- I2X and then average over x without altering the value of 
the sum 

^2 E ^Ex"(^^2/i)x:((Q2/2)x(Qr'^r')2;)- m 

To evaluate the group average over x we now use the identity X^^gc Xo(«5^5) = 'j~Xa{o-b^^) 
(see Eq. (10) with Cq = 1 for real representations). This gives 

7"j7^ E XaiPl2h)X*a{Q2l2hQl) = E Xa{9p)X*ai9q') ■ (61) 

Hence by performing an antiparallel encounter exchange en route to changing B L) D to AU C we 
have expressed k^™''* by ^ times the coefficient for a structure where q is replaced by q', i.e. 

= {Xa{9pi)Xa{gp2) ■ ■ ■ Xa{9p)X*a{9qi)X*ai9q2) ■ ■ ■ X*a{9q')) ■ (62) 
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If after the reconnection p and q' are identical then we have gp = gqi allowing us to reduce the 
number of variables further since {Xa[gp)x*a{9p)) — 1- 

We must now show that performing all other types of reconnection steps leads to the same 
result, allowing k^'"'" = Sa^^ to be infered by recursion. In one alternative situation B U D 
will contain l^^ and I2 in two separate orbits indicating that the relevant antiparallel encounter 
stretches also belong to different orbits. However since the equality X%{9~^) = Xaig) = Xaid) (foi' 
real representations) allows us to invert any of these periodic orbits without fear of altering the 
group sum we can convert this into a situation with parallel encounter stretches and non-inverted 
link elements li and I2 as treated below. 



5.1.2 Parallel encounter permutations 

We now deal with the case where BUD contains the stretches h and I2 without inversion; the case 
of and I2 appearing can be reduced to this situation by reverting all orbits in time. There are 
two scenarios: Either li and I2 and thus the relevant encounter stretches in B L) D belong to the 
same periodic orbit {q say) or they belong to different orbits {qi and q2 say), illustrated in Figs. [7] 
and [8] respectively. In the former case performing the encounter exchange will split the orbit into 
two separate orbits (denoted by q[ and q'2), see Fig. ff\ whereas in the latter case the two orbits will 
merge into a single orbit (denoted by q'), see Fig. |8| 



e, 

, > 



Figure 7: (a) The permutation of parallel encounter stretches within a period orbit q leads to (b) 
a decomposition of the orbit into two orbits q'^ and (^2- After the permutation the orbit q has an 
encounter stretch that is almost identical to a stretch of the orbit p. 

We first consider the case where li and I2 are contained in the same orbit q £ B U D. In this 
case the group elements of p and q can be written in the form 

gp = Phh and Qq = {Q2){hQih)- (63) 

As ^1 and I2 follow each other in p the orbit q must contain almost parallel encounter stretches 
following the link with group element li and preceding the link with group element l2- The positions 



of these encounter stretches are indicated by the brackets in the factorisation of gq in (63). Now 
switching connections between the encounter stretches leads to a decomposition of q into two orbits 
(^i and ^2 with group elements 

9q[ = hQih and 5^/ = Q2. (64) 

As the links with group elements li and I2 follow each other in q'^^ just as in p we have thus performed 
one of the reconnection steps leading from B [J D to A[J C . As before we now relate the coefficient 
^struct £qj, q^j, original structure to the one for the structure with q replaced by q'l and (^2- The 
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coefficient k^''"''' involves tfie average 

^ Xa{Phh)X*a{Q2hQll2) 



\G\' 



(65) 



We now make tlie same substitutions h ^ x ^li and I2 — )• hx as in the previous case to set up an 
average over x of tlie form 



j?^ T?T\YXa{Phh)X*aiQ2X HiQihx) 



(66) 



Using the identity X^geG Xa{agbg ^) = -^Xa{a)xa{b) (see Eq. (12)) as well as the invariance of 
the characters under cyclic permutation of the group elements we men obtain 



1 1 



Xa{Phh)x:AhQih)xl{Q2 



1 1 



Yj Xa{9p)X*a{9q[) Xa \ 



(67) 



If we insert this result into k!^''"''* we obtain once again a factor of — times the contribution of a 
structure where one reconnection step has been performed to make B U D more similar to ^ U C, 
i.e. 



1 



Xaigpi)Xa{gp2) ■ ■ ■ Xa{gp)X*a{9qi)X*ai9q2) ■ ■ ■ X*ai9q[)X*ai9< 



^12' 



(68) 



We should mention that there can be a special case in which the orbit p has only one encounter 
stretch and one link, in which case the two links h and I2 would be identical. In this situation 



(63) is replaced hy gp = I and gg = Ql, and we arrive at the same final result (68) if we make the 
substitution / — t- xlx^^. 



pe A(jC 



9i 



Figure 8: (a) Exchanging two parallel encounters belonging to separate encounter stretches leads 
to (b) the recombination into one orbit of the combined length. 

Now let us assume that the links li and I2 belong to two different orbits {qi and ^2) in B U D 
as depicted in Fig. [8} In this case the group elements have the form 

gp = Phh, gqi=hQi and gq2 = Q2h- (69) 

By the same arguments as above there are almost parallel encounter stretches following the link 
element associated to li and preceding the link element associated to Z2. If we switch connections 
between these stretches the orbits qi and q2 merge into a single orbit with the group element 

gq' = {Q2l2){hQl). (70) 
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Now the coefficient k^''"''* involves the average 



^ Xc.{Phh)xl{Qlh)xl{hQ2) = V^, Xa{Pl2h)X*a{hQlX-')x*ail2xQ2), (71) 

' ' h,h&G ' ' x,h, 



where we have used the invariance of the characters under cyclic permutation of group elements 
and the same substitutions li — )• x~^li and I2 — )■ hx as in the previous two cases. We then use the 
identity tL Y.g^GXa{ag)Xa{hg'^) = j^Xa{ab) (see Eq. 11) and sum over x to get 



1 1 



1(712 X] Xa{Pl2h)X*aihQlQ2h) = 77^ Xa{gp)x*ai9q 
" ' ' ll,l2&G " ' ' h,l2&G 



(72) 



^struct 



is given by 



Thus after exchanging encounter stretches we again obtain the result that kJ, 
times the contribution of a structure with one reconnection step removed. 

We thus see that performing each of the L — V steps needed to transform the set of periodic 
orbits in B U D into AU C produces a factor of s~^. When one of the orbits in B U D has been 
turned into one already appearing in A U C this orbit can be dropped from the sum in k^''"''* due 

1. After all L — V reconnections have been performed no more orbits 
for this trivial structure is just 1. Altogether this implies the desired 



to iG\^9eGXaig)xU9) 



3 lei 
result 



are left and the value k^'"'" 



SO Eq. (50) simply becomes 



.(1) 



■,off 



struct 



(1) 

off, struct 



(73) 



(74) 



This sum was already shown to yield GOE behavior in [22j. Hence the correlation function Ca{e) 
for subspectra associated to real representations in time reversal invariant systems is faithful to the 
GOE prediction, obtained from Eq. ([s]) after dropping the restriction to the real part. 



5.2 Complex Irreducible Representations 

If the representation is complex then the difference between parallel and anti-parallel exchanges be- 



comes greatly significant. For the parallel exchanges in section 5.1.2 nothing changes in comparison 



to the real case, since the group relations (11) and (12) hold regardless of the type of represen- 



tation. However, for the antiparallel exchanges in section 5.1.1 the group relation (10) implies 



a zero contribution when the representation is complex. The underlying reason for this rests on 
the inequivalence between a complex representation and its complex conjugated counterpart: An 
antiparallel exchange implies the existence of a link, say la, which is traversed in opposite directions 
in j4 U C and B U D, corresponding to a time-reversal. If we express the characters as traces of 
the matrices M^'^^g) and use the representation property M^°'\g2gi) = M^°'\g2) M^°'\gi) then it 
becomes clear that the average over this one link involves a calculation of the type 



|G| ^ 



Ik i^a) 



0, (75) 



where we used in the first step the unitarity of the representation, in the second step /3 denotes 
the complex conjugate representation of a, and in the third step we used the group orthogonality 
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relation ([s]). The contributions of antiparallel exchanges that in the real case could be converted 
to parallel exchanges by reverting orbits in time vanish for the same reason. Hence only those 
structures contribute where the periodic orbits in BDD follow all links and encounters in the same 



direction as the periodic orbits in ^4 U C. Thus we obtain a sum of the form (74) but restricted to 
structures which do not require time-reversal invariance. This sum was already shown in [22] to 
yield the GUE behavior. 

It is interesting that the structures involving time reversal are dropped not because they are 
non-existent but because over many pseudo-orbits their average contribution is zero. 

5.3 Non Time-Reversal Invariant Systems 

For systems without time reversal invariance the structures relying on time reversal are obviously 



excluded, this time due to the dynamics of the system. We obtain again the result (50) for the 



off-diagonal contributions, but the sum is now restricted to structures that do not require time- 



reversal invariance. The calculations in section 5.1.2 on parallel encounter permutations carry over 
unchanged to show that k^^^"^^ = s^^^. As mentioned before, these calculations do not depend 
on the type of the representation, and in this way one obtains GUE behaviour for both real and 
complex representations. 



6 Full Correlation Function 

To describe the correlations in the full quantum system we need to understand not only correlations 
inside each subspectrum a but also correlations between different subspectra a and /3. These cross- 
correlations were also considered in [23] where the leading terms in ^ as well as the first subleading 
non-oscillatory term were determined. 

When defining the cross-correlation function between the subspectra a and /3 one has to take 
into account that the average level densities for the two subspectra can be different, and one could 
scale the energy according to the average level densities for either of the subspectra or for the full 
system. The corresponding dimensionless parameters e, e*-"\ and e^^-* are related by | = = 

For the moment we will use e*-"^ and e*-^^ in parallel and define the cross-correlation function by 

Ca^ie^'^^e^^^) = ( -trP„ 1 -trP, i ) - J . (76) 

We consider first the case when a and (3 are complex conjugate representations. 



6.1 Scenario 1: a and (5 mutually complex conjugate 

The subspectra associated to two mutually complex conjugate representations are identical in time- 
reversal invariant systems, as was discussed in section [2] As a consequence, their cross-correlation 



function (76), with e = e*-"^ = e^^-*, must coincide with the correlation function Ca(e) for complex 
representations, i.e., with the GUE result ([2]). This can also be seen from semiclassics if we use 
that two mutually complex conjugate representations a,/3 satisfy Xpio) = Xaig) = Xa{g~^)- Hence 
the character associated to an orbit in the representation /5 coincides with the character associated 
to the time-reversed orbit in a. This means that whenever the representation /3 is considered we 
have to work with time-reversed orbits instead. As a consequence the result from the diagonal 
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approximation (37) is the same as before because it is unchanged if partner orbits are replaced 



by their time-reversed versions. Concerning the off-diagonal contributions, the relevant structures 
coincide with the structures used before apart from time reversal of all orbits associated to /3. 
These structures are in a one-to-one relation to those considered previously and give the same 
contributions. 

For systems without time-reversal invariance mutually complex eigenfunctions are not required 
to possess the same energy and so the two-fold degeneracy occurring between mutually-complex 
representations disappears. From a semiclassical perspective, this result arises because there are no 
time-reversed orbits in our system, and the contributions from parallel encounters vanish because 



of similar arguments as in ( 75 ) . 



6.2 Scenario 2: a and /3 not mutually complex conjugate 



We now consider a and /? that are neither identical nor mutually complex conjugate. We will show 
that the corresponding subspectra are uncorrelated in the semiclassical limit. However this does 
not rule out correlations outside the semiclassical limit as observed in [30] . The generating function 



that corresponds to the cross-correlation function (76) is defined as 



.(a) (/3) (a) (/3) 



A. 



E 




(77) 



and gives the cross correlation function through 



(ID 



(78) 



where 



implies 



-B 



The part of the generating function 



responsible for non-oscillatory contributions is given by 



^(1) 




G^-* G^"** G^''* ( — 1 ) """"^ 



,i(S^-SB+Sc-So)/ri i(TA.if^+Tce^"')/Tf,-' ilTB^'^^+To^^O/rif' 



(79) 



generalizing Eq. ( 29 ) . Again the action difference becomes small if the orbits in B U D coincide 



with those in ^ U G apart from their connections inside encounters. However the orbits in A U G 
contribute with their character in the representation a whereas the orbits in BUD contribute with 
the complex conjugate of their character in the representation /3. Hence all contributions vanish 
due to PI X^geG ~ ^' particular this relation comes into play in Eq. (34) in the 

diagonal approximation and when dealing with orbits that have become identical after the recursion 



steps in our calculation of k^''"'''. Thus the only contribution to (79) arises from the empty set, 
where = ub = nc = n^) = 0. This summand trivially gives 1 and together with the oscillatory 
prefactor we have 



Z 



(1) 

a/3 



■1 {") 



(80) 



21 



(2) 

The second part Z^^ of the generating function has to be considered separately, since the inclu- 
sion of different representations mean it can no longer be obtained directly from Z^^ through Eq. 



(26). Instead the complex conjugation of both spectral determinants in the numerator, following 



Riemann-Siegel resummation, yields 

\a,b,c,d 

The action difference now becomes small if the orbits in i? U C are related to those in Au D. 
However if a ^ (3 and there is (part of) an orbit in A which is correlated to (part of) an orbit in B 
then the group orthogonality relation for different representations will instill a zero contribution. 
The same process occurs for those (parts of) orbits in C which are correlated to (parts of) orbits 

(2) 

in D. The generating function may therefore be split into two uncorrelated factors confined 
by the two irreducible representations 



= ei(4"'+4"'*)/2 /^Gj)Gj.")*(-l)"Ce^(^-^-^c)/'"'e*(^^^A^-T^c4"^*)/Tif) \ 
\a,c I 

X e*(^s'+4^'*)/2 / ^G^)*G(J)(_i)nDe-*(^s~^^)/V(^s^s'-^D4^>*)/<'\ . (82) 
\b,d I 

f2) ('2) 

This can once again be separated into the form Z\^1^{1 + Z^^ ) where the diagonal part is given by 



ap,diag 



\ a 

'exp 



X e^(4''+4''*)/2exp -4-y|Fj2e*7^44''-4''*)/TifM . (33) 



If we let X = ^(e^^ — e^T*^*), then Rex > 0, and after invoking Hannay and Ozorio de Almeida's 



sum rule [9] the sum over a in the first line of ( 83 ) turns into 



^|^^|2,«T/T("' ~ r^e^-T/T'-^ = r ^ ^e-- (84) 

and hence diverges in the semiclassical limit due to — t- c«. The same applies for the sum in 

(2) 

the second line and altogether we have Z^^ ^.^^ — )• 0. The off-diagonal part is given by 



.(2) _ ^ (- i)"gs^-^Kr" 1 

2VV\{-1)L (,,W_,^W*)L-y 



7(2) _ V- 



-c 



( -l\nj:) -L—V struct ^ 



struct ^ ' \'''-B ^ D 
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which remains finite even when = e[?^ or = e'^ in the denominator. Thus Z^^^ can indeed 
be dropped leaving only 



Z^p = Z^^l = e<'A +-B )/2 



and thus combined with relation (78) 



(86) 



(87) 



(li) 



Hence in the semiclassical limit there are no correlations between subspectra associated to different 
representations unless these representations are mutually complex conjugate. 



6.3 Final result 

The correlation function of the full system is given by a linear combination of all correlation 
functions inside and between subspectra. Taking into account the scaling factors from p = s^Pa 



and e^"-* 



\G\ 



e we obtain 



SaPaSi3Pl3 



Ca/3(e 



) = E 



For time-reversal invariant systems we have shown that 



Ca/3(e 



'CGOE(e^"^) ifa = /3real 

CGUE(e^"^) if a = f3 complex or a, /3 mutually complex 
otherwise 



yielding 



C{e) 



1 

w 



a real 



GOE \ iQi^ j + E ^-^a^GUE ( TT^' 
a complex 



For systems without time-reversal invariance we have shown that 

/CGUE(e(")) if a = 13 



C.,(6(-),6(^)) 



otherwise 



and thus 



1 



C(e) — E ^Q^GUE f VTZVe 



\G\ 



(88) 



(89) 



(90) 



(91) 



(92) 



This is our final result: It shows that a time-reversal invariant quantum chaotic system has a 
correlation function in the semiclassical limit that is comprised of sums of GOE and GUE correlation 
functions corresponding to real and complex subspaces and weighted by appropriate degeneracy 
factors (provided that the fundamental domain does not possess any hidden symmetries leading to 
deviations from RMT as in arithmetical billiards [41j). In contrast, a non-time reversal invariant 
system has no additional degeneracy factor due to complex conjugate representations and every 
subspace has a GUE correlation function. 
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7 Conclusions 



For chaotic systems with a discrete spatial symmetry we have shown that periodic-orbit theory can 
account for correlations of levels inside each of the subspectra that are associated with the different 
irreducible representations a of the symmetry group as well as for cross-correlations between these 
subspectra. This was achieved by considering the semiclassical realisation of the symmetry-reduced 
spectral determinant Aa{E) and implementing the Riemann-Siegel lookalike formula, from which 
the exact asymptotic RMT expansions for the complex correlation functions were successfully 
reproduced. The results are summarised in Table [T] For real representations in time-reversal 
invariant systems we obtained GOE behaviour and otherwise GUE behaviour. 





no T inv. 


r2 = i 


Q complex 


GUE 


GUE 


a real 


GUE 


GOE 



Table 1 : Table showing how the random matrix ensembles for subspectra depend on representation 
type as well as time reversal properties of the system. 

Throughout this paper we have neglected the third type of irreducible representations, known as 
pseudo-real (with complex representation matrix M^"^ but real trace Xa)- These arise for symme- 
tries which are more complex than the standard crystallographic point groups, with the simplest 
being the quaternion group Q8. However, from a semiclassical analysis this case is interesting; 
showing, contrary to earlier predictions [4j, to possess a GSE distribution rather than GOE [30]. 
We have also neglected the case of time-reversal operators squaring to minus unity r2 = -1 as IS 
the case in systems with half-integer spin. For these systems we expect GUE for complex repre- 
sentations, GSE for real representations, and GOE for pseudo-real representations. Including these 
cases would extend Table [T] to a 3 x 3 table. 

The techniques in this paper could also be used to analyse the effects of false time reversal 
symmetry [12] ) continuous symmetries [53] , symmetry breaking [441 HSj BH] , and arithmetical sym- 
metries such as the eigenvalues of the Laplace-Beltrami operator on the domain of the modular 
group [41j. 
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